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Equations  have  been  derived  for  determining  the  dynamic  response  of  the 
combination  of  a  linear  elastic  structure  and  its  subsystems.  The 
method  is  based  on  measured  mobility  matrices  of  the  primary  structure 
and  the  subsystem  independently.  Mathematical  relationships  were  for¬ 
mulated  for  the  main  system  and  subsystem  interface. 

The  mathematical  model  established  provides  for  a  wide  range  of  cross 
coupling  effects  to  simulate  diverse  subsystems.  Specifically,  the 
types  of  subsystems  considered  were  a  spring-mass  system  connected  at 
one  point,  a  rigid  inertial  mass  elastically  connected  at  two  points, 
and  a  beam  elastically  attached  at  three  or  more  points.  The  spring- 
mass  subsystem  attached  at  one  point  is  illustrative  of  a  single  vibra¬ 
tion  absorber  or  a  load  suspended  from  the  helicopter.  The  rigid 
inertial  mass  subsystem  with  two-point  attachment  is  typical  of  a 
munitions  store  or  fuel  storage  tank.  The  beam-type  subsystem  connected! 
at  three  or  more  points  is  representative  of  a  suspended  weapon  or 
possibly  an  auxiliary  engine. 

A  digital  computer  program  vis  generated  for  the  IBM  Model  360/40 
computer  using  FORTRAN  IV  language  to  numerically  test  the  afore¬ 
mentioned  theory.  Computer  experiments  were  conducted  to  test  tne 
sensitivity  of  the  theory  to  measurement  error  in  the  simulated  rest 
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FOREWORD 


The  work  presented  in  this  report  was  performed  by  Kaman 
Aerospace  Corporation  under  Contract  DAAJ02-70-C-0012 
(Task  1F162204AA4301)  for  the  Eustis  Directorate,  U.  S.  Army 
Air  Mobility  Research  and  Development  Laboratory,  Fort 
Eustis,  Virginia.  The  program  was  implemented  under  the 
technical  direction  of  Mr.  Joseph  H.  McGarvey  of  the 
Reliability  and  Maintainability  Division* and  Mr.  Arthur  J. 
Gustafson  of  the  Structures  Division  .**  The  report  is  pre¬ 
sented  in  four  volume  s ,  each  describing  a  separate  phase  of 
the  basic  theory  of  structural  dynamic  testing  using  impedance 
techniques . 

Volume  I  presents  the  results  of  an  analytical  and  numerical 
investigation  of  the  practicality  of  system  identification 
using  fewer  measurement  points  than  there  are  degrees  of 
freedom.  The  parameters  in  Lagrange's  equations  of  motion, 
mass,  stiffness,  and  damping  for  a  mathematical  model  having 
fewer  degrees  of  freedom  than  the  linear  elastic  structure  it 
represents  may  be  determined  directly  from  measured  mobility 
data.  Volume  II  describes  the  method  of  system  identification 
wherein  the  necessary  impedance  data  are  experimentally  deter¬ 
mined  by  applying  a  force  excitation  at  a  single  point  on 
the  structure.  Volume  III  presents  a  method  of  determining 
the  free-body  dynamic  responses  from  data  obtained  on  a  con¬ 
strained  structure.  Volume  1^  describes  a  method  of 
obtaining  the  equations  for  the  combination  of  measured 
mobility  matrices  of  a  helicopter  and  its  subsystems.  The 
response  of  the  combination  of  a  helicopter  and  its  sub¬ 
systems  is  determined  from  data  based  on  the  experimental 
results  of  the  main  system  and  subsystems  separately. 


**Division  name  changed  to  Military  Operations  Technology  Division 
’^Division  name  changed  to  Technology  Applications  Division. 
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INTRODUCTION 


The  success  of  a  helicopter  structural  design  is  highly 
dependent  on  the  ability  to  predict  and  control  the  dynamic 
response  of  the  fuselage  and  appended  components.  An  effec¬ 
tive  dynamic  analysis  of  complex  systems  should  yield  the 
response  of  the  primary  system  and  its  associated  subsystems. 
It  is  extremely  desirable  to  analytically  predict  the  com¬ 
plete  response  of  a  linear  elastic  structure  due  to  the 
addition  or  alteration  of  particular  components. 

This  report  describes  a  method  whereby  the  dynamic  response 
of  the  entire  system  can  be  determined  from  knowledge 
of  the  response  of  the  main  system  and  the  subsystem  sep¬ 
arately.  The  formulation  is  predicated  on  the  theory  of 
structural  dynamic  testing  using  impedance  technigues.  The 
analysis  requires  measured  mobility  matrices  for  the  basic 
structure  alone  and  free  mobility  matrices  for  the  attached 
component.  Therefore,  once  the  mobility  matrices  for  the 
basic  system  are  measured,  they  can  be  continually  used  in 
conjunction  with  the  measured  free  mobilities  for  the  various 
components  connected  to  the  main  structure . 

It  is  anticipated,  in  practice,  that  the  mobility  matrix  for 
the  basic  system  would  be  obtained  by  the  method  of  Volume  II 
of  this  report  and  that  the  free  mobilities  for  the  components 
would  be  obtained  by  the  method  of  Volume  III. 

Specifically,  in  the  present  report  the  primary  system  was 
a  20-degree-of-freedom  representation  of  an  actual  helicopter, 
and  three  types  of  subsystems  were  considered.  The  sub¬ 
systems  studied  included  a  spring-mass  system  connected  at 
a  single  point,  a  rigid  inertial  mass  elastically  attached 
at  two  points  ,  and  a  beam  elastically  connected  at  three  or 
more  points.  The  method  employed  simulated  test  data  to 
represent  the  required  experimental  mobility  data, and  measure¬ 
ment  errors  were  introduced  to  test  the  sensitivity  of  the 
theory  to  error. 


THEORY 


Consider  a  finite  degree  of  freedom  simulation  of  an  actual 
helicopter.  The  impedance  matrix  for  the  system  can  be 
expressed  in  terms  of  mobilities 
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The  impedance  of  a  subsystem  to  be  attached  to  the  primary 
system  can  also  be  expressed  in  terms  of  mobilities 


(2) 


The  mobility  of  the  combined  system  is  defined  as  the  inverse 
of  the  impedance 
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The  product  of  the  impedance  matrix  and  the  mobility  matrix 
is  the  unit  matrix 


[zi  [y]  =  ([y]"1  +  t^r1)  [YI  “  II] 


(4) 


Multiplying  both  sides  of  Equation  (4)  by 
[Y]  yields 

[Y]  =  ([I]  +  [YHYg]"1)  "1[Y’ 


[Y]  and  solving  for 
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Substituting  the  actual  matrices  into  the  matrix  Equation  (5), 
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Performing  the  indicated  operations  within  the  matrix  inverse 
results  in 
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To  facilitate  solution  of  the  matrix  Equation  (7) ,  the  in¬ 
verse  of  the  matrix  on  the  right  side  of  the  equation 
must  be  evaluated.  Let 
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Solving  for  the  elements  of  the  X  matrix  yields 
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Substituting  the  X  matrix  as  expressed  in  Equation  (10)  into 
Equation  (7)  and  expanding  gives 


(12) 

Equation  (12)  yields  the  response  of  the  complete  system 
for  force  excitation  at  each  position  on  the  main  system 
including  the  interface  points.  Since  the  force  excitation 
on  a  helicopter  i_.  usually  applied  at  one  particular  point, 
Equation  (12)  can  be  reduced,  yielding  the  complete  struc¬ 
tural  response  for  forcing  at  a  single  point  on  the  structure. 
.Thus,  Equation  (12)  reduces  to  a  column  of  mobilities  for 
force  excitation  applied  at  position  j  of  the  structure 


(13) 


where  YftA  represents  the  dynamic  response  of  each  position 
on  the  ■'primary  system  excluding  the  interface  points  and 
[YBAj describes  the  response  of  the  attachment  points. 

It  is  possible  to  obtain  the  dynamic  response  for  the  complete 
system  utilizing  the  modal  matrix  of  the  points  of  interest  on 
the  primary  structure  exclusive  of  attachment  points ,  the 
modal  matrix  of  the  interface  points  and  the  complex  modal 
mobility  of  the  primary  system  and  the  free  mobility  of  the 
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appended  subsystem.  These  parameters  can  be  obtained 
employing  the  techniques  described  in  References  1  and  2. 
Following  the  method  of  Reference  2,  the  mobility  of  the 
main  system  exclusive  of  the  subsystem  connection  points 
is  given  by 


IYAB]  -  [\] 


(14) 


where  the  number  of  rows  of  the  matri  corresponds  to 
the  number  of  points  of  interest  on  the  nain  system,  ex¬ 
clusive  of  the  subsystem  attachment  points ,  and  the  number 
of  columns  corresponds  to  the  number  of  modes  and  is  equal 
to  the  total  number  of  points  of  interest  on  the  primary 
system  including  the  interface  points.  Matrix  £  |  has  the 

same  number  of  columns  as  the  number  of  connection  points 
between  the  two  systems  and  the  identical  number  of  columns 
as  matrix  £  1.  The  mobility  of  the  subsystem  attachment 
point  is  L  a] 


[YBB]  =  [*BirY*JIV 


(15) 


If  Equations  (14)  and  (1.5)  are  substituted  in  Equation  (13), 
the  result  is 


ERROR  ANALYSIS 


(16) 


Measurements  of  the  complex  mobilities  will  be  subject  to 
experimental  errors  of  various  types  ,  including  errors  in 
equipment  calibration,  errors  resulting  from  equipment 
incompatibility,  errors  due  to  extraneous  signals  and  errors 
due  to  random  noise. 

Generally,  all  errors  can  be  classified  as  either  random 
or  bias.  The  random  errors  are  equally  likely  to  be 
positive  or  negative,  whereas  the  bias  errors  are 
systematic  and  in  one  direction  only.  In  the  present  study, 
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both  types  of  measurement  error  have  been  included.  The 
simulated  experimental  data  were  polluted  with  measurement 
errors  of  +5  percent  random  and  5  percent  bias  or  both  the 
real  and  imaginary  components  of  displacement  mobility. 

As  indicated  in  Reference  3,  there  is  no  definitive  prob¬ 
ability  distribution  for  errors  oi:  each  type  in  impedance 
testing  practice.  In  the  present  analysis,  a  random  number 
generator  was  utilized  with  a  resultant  uniform  distribution 
o2  random  error.  The  rectangular  distribution  of  accidental 
type  error  between  the  selected  limits  is  very  conservative 
compared  to  the  usual  definition  of  the  limits  at  hree 
standard  deviations  from  the  mean  of  a  normal  distribution. 
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COMPUTER  SIMULATION  RESULTS 


A  computer  study  to  determine  the  response  of  the  combination 
of  a  helicopter  and  its  subsystems  based  on  the  simulated 
test  results  of  the  individual  system  and  subsystems  can  be 
extremely  useful  in  the  development  cycle  of  a  helicopter. 

In  the  present  analysis,  a  mathematical  model  was  established 
to  provide  for  a  wide  range  of  cross*- coupling  effects  to 
simulate  diverse  subsystems.  The  helicopter,  or  main  system, 
was  represented  by  a  20-degree-of-freedom  mathematical  model. 
Table  I  presents  a  lumped  mass  description  of  the  afore¬ 
mentioned  specimen  which  was  used  to  generate  the  simulated 
experimental  data.  Three  types  of  subsystems  were  incor¬ 
porated  in  the  study,  represented  as  a  spring-mass  system 
elastically  connected  at  one  point,  a  rigid  inertial  mass 
elastically  connected  at  two  points  and  a  beam  elastically 
connected  at  three  or  more  points.  Figure  1  shows  the  afore¬ 
mentioned  subsystems. 

Figures  2  and  3  present  the  real  and  imaginary  displacement 
mobility  frequency  response  for  the  main  system-subsystem 
interface  point  for  a  spring-mass  subsystem.  The  force 
excitation  was  applied  at  Station  6,  the  hub  ptation,  and 
the  response  was  measured  at  Station  3,  the  general  area  of 
the  pilot  seat,  and  coincident  with  the  subsystem  attachment 
point.  Data  are  presented  for  conditions  of  zero  experimental 
error  and  for  simulated  experimental  displacement  mobility 
data  recorded  with  a  random  error  of  +5  percent  and  a  bias 
error  of  +5  percent.  For  the  cases  involving  error,  the 
random  displacement  mobility  error  was  computed  using  a 
uniformly  distributed  probability  density  function.  This 
error  was  applied  to  both  the  real  and  imaginary  components 
of  the  main  and  subsystem  displacement  mobility  data.  As 
can  be  observed  from  Figures  2  and  3,  the  method  is  extremely 
insensitive  to  the  measurement  error  as  applied  herein. 

Figures  4  and  5  show  the  same  type  data  as  given  in  the 
previous  figures  except  that  the  subsystem  investigated  was 
a  rigid  inertial  mass  elastical.ly  attached  to  the  main 
structure  at  two  points.  The  interface  points  in  this 
situation  are  located  at  Stations  1  and  2  of  the  main  system. 
Again,  the  frequency  response  of  the  displacement  mobility, 
both  real  and  imaginary,  is  effectively  invariant  with 
error  for  the  error  level  incorporated.  Figures  6  and  7 
present  the  same  type  results  for  the  beam  subsystem  elas¬ 
tically  connected  at  three  points.  The  attachment  points 
in  this  instance  are  at  Stations  1,  2,  and  3  of  the  main 
system.  The  data  substantiate  the  previous  observations 
of  the  relative  insensitivity  of  the  method  to  error. 
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SPRING-MASS  SUBSYSTEM 


RIGID  INERTIAL  MASS  SUBSYSTEM 


K1  =  K2  »  K3  =  Kn  =  4000  LB/IN. 
M  =  2.15  LB-SEC2/IN. 


| _  -  120" - 

BEAM  SUBSYSTEM 

Figure  1.  Subsystem  Representation. 


9 


N 


DISPLACEMENT  MOBILITY,  REAL  r IN. /LB- FORCE) 


FREQUENCY  -  HERTZ 

Figure  2.  Real  Displacement  Mobility  Frequency  Response; 
Combination  of  Main  System  and  Spring-Mass 
Subsystem. 
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cjid  Inertial  Mass  Subsystem. 


SYMBOL  SEED  ERROR 


14 


The  effect  of  a  spring-mass  subsystem  on  the  response  of  the 
20-degree-of-f reedom  mathematical  representation  of  the  heli¬ 
copter  is  shown  in  Figure  8.  A  sinusoidal  force  excitation 
was  applied  at  the  hub  station,  and  the  real  displacement 
mobility  is  given  for  the  basic  system  alone  and  for  the 
main  system  with  subsystem  attached  at  Station  3,  the  pilot 
seat  location.  The  forcing  frequency  was  9.2  Hertz,  which  is 
the  calculated  natural  frequency  of  the  spring-mass  subsystem 
alone.  As  would  be  expected,  the  response  of  the  total 
system  is  characterized  by  a  nodal  point  at  the  interface 
station.  The  effect  of  error  is  also  indicated  on  the  figure. 
The  dispersion  of  the  results  using  the  various  random 
error  seeds  is  within  an  acceptable  level.  Figure  9 
presents  similar  data  as  Figure  8  except  that  the  forcing 
frequency  is  20  Hertz.  The  disturbing  frequency  is 
separated  enough  from  the  subsystem  resonant  frequency  that 
there  is  essentially  no  difference  in  the  system  response 
with  or  without  subsystem  attached.  Figures  10  and  11 
illustrate  the  influence  of  the  rigid  inertial  mass  subsystem 
on  the  response  of  the  combination  of  the  basic  helicopter 
representation  and  the  appended  subsystem.  The  forcing  fre¬ 
quencies  used  represent  the  approximate  natural  frequencies 
of  the  subsystem  alone.  The  effect  of  error  can  be  observed 
from  the  figures,  and  the  insensitivity  of  the  analysis  to 
measurement  error  is  again  visible.  The  response  of  the 
configuration,  consisting  of  the  main  system  and  a  beam  elas¬ 
tically  connected  at  three  points  on  the  main  structure,  is 
represented  in  Figure  12  for  an  excitation  frequency  of  10 
Hertz  and  in  Figure  13  for  a  forcing  frequency  of  30  Hertz. 

The  figures  also  yield  the  effect  of  measurement  error  on 
the  system  displacement  mobility  response. 

The  results  shown  represent  a  small  number  of  the  computer 
simulation  experiments  actually  conducted.  For  each  of  the 
subsystem  configurations  considered,  the  location  of  the 
attachment  points  along  the  main  system  could  be  varied  :s 
long  as  compatibility  between  the  systems  was  maintained. 

The  flexibility  of  the  method  presented  in  this  report  in 
determining  the  response  of  the  combination  of  a  helicopter 
and  its  subsystems  based  on  test  results  of  the  individual 
system  and  subsystems  provides  a  means  of  modifying  subsystem 
characteristics,  interface  locations  and  connection  effects. 

A  wide  range  of  cross-coupling  effects  to  simulate  diverse 
subsystems  can  be  analyzed  based  on  the  measured  mobilities 
of  the  individual  subsystems,  since  the  mobility  data  of  the 
main  system,  once  measured,  remain  constant. 

A  digital  computer  program  listing  in  FORTRAN  IV  lancruage 
and  a  description  of  the  program  input  cards  are  given  in 
the  appendix. 
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Figure  10.  20-Degree-of-Freedom  Model;  Real  Displace: 

Effect  of  Rigid  Inertial  Mass  Subsystem. 
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Figure  11.  20-Degree-c f-Freedom  Model;  Real  Displacement  Mobility;  With  and 
Without  Rigid  Inertial  Mass  Subsystem.  Frequency  =  28  Hertz. 


CONCLUSIONS 


1.  The  response  of  the  combination  of  a  helicopter  and 
its  subsystems  can  be  determined  based  on  the  test 
results  of  the  individual  system  and  subsystems. 

2.  The  method  is  insensitive  to  measurement  error  using 
simulated  test  data  subjected  to  errors  that  are 
within  the  state  of  the  measurement  art. 

3.  The  method  can  be  used  to  study  a  wide  range  of  cross 
coupling  effects  to  simulate  diverse  subsystems. 

4.  At  a  specific  excitation  frequency,  once  the  mobility 
data  for  the  main  system  are  measured,  they  remain 
constant;  thus,  only  the  measured  mobilities  of  the 
individual  subsystems  need  be  considered. 

5.  The  method  provides  an  expedient  means  of  modifying 
the  subsystems  appended  to  a  helicopter  at  a  develop¬ 
ment  stage  of  the  system. 

6.  The  method  is  amenable  to  experimental  implementation 
and  is  numerically  sound. 
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APPENDIX 

COMPUTER  PROGRAM  DESCRIPTION 


A  digital  program  was  prepared  for  determining  the  response 
of  the  combination  of  a  helicopter  and  its  subsystems  based 
on  the  test  results  of  the  individual  system  and  subsystems. 
The  program  was  written  for  the  IBM  360/40  disk  operating 
system  using  FORTRAN  IV  language.  A  flow  chart  depicting 
the  program  logic  is  shown  in  Figure  14 .  A  description  of 
the  input  cards  and  a  program  source  listing  are  presented 
in  this  appendix. 
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DESCRIPTION  OF  INPUT  CARDS 


Note :  All  integer  variables  must  be  right  justified  with 
no  decimal  point. 

Tape,  Cat d  Reader  and  Printer  Assignments 
1  Card  Reader 
3  Pointer 

11  Contains  free  displacement  mobility  matrices  for 
subsystem  connection  points. 

12  Contains  displacement  mobility  matrices  for  primary 
system. 

All  input  data  must  be  in  the  following  units: 

2 

Mass  -  lb-sec  /in. 

Stiffness  -  lb/in. 

Frequencies  -  Hz 
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PROGRAM  COMSYSA 
COMPONENT  SYNTHESIS 


Card 

1  Columns  1-10 

M 

Number  of  interface  points 
(Number  of  attachment  points 
between  subsystem  and  primary 
system  (FORMAT  110). 

MS 

Interface  points,  station 
numbers.  10  columns  per 
value  (Maximum  10  values). 
(FORMAT  8110) . 

NC0L 

Primary  system  station 
number  at  which  forcing 
function  is  applied  (FORMAT 
110)  . 

Card 

2  1-10 

PCTR 

Random  error  applied  to  real 
mobility.  Uniform  between 
-  and  +  PCTR*  ELEMENT 
displacement. 

11-20 

PCTBR 

Bias  error  applied  to  real 
mobility.  PCTBR*  ELEMENT 
displacement . 

21-30 

PCTI 

Same  as  PCTR  except  applied 
to  imaginary  displacement 
mobility. 

31-40 

PCTBI 

Same  as  PCTBR  except  applied 
to  imaginary  displacement 
mobility. 

41-50 

IZ 

Random  error  seed,  used  in 
generation  of  error. 

Card (s ) 

3  1-10 

KEEP 

Stations  to  be  used  in 
primary  system,  co  columns 

per  value,  8  values  per  card 
(FORMAT  8110) . 
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Jg^SiaiSsSaB-i. 


c 

c 


c 


COMPONENT  synthesis 

-NP 

1 

MNP 

2 

DIMENSION  MSI  201. HMTI  .HEA0I201  ,»U1  iUu  1 ,  YR  ( 20,  211  .VI120.21I. 

MNP 

3 

AVABR (2U,2i) «YABI(20,2l),  YCCRI2j,«.1J , YcCI  1 20,2 1 1 ,V  AAAI 20,2 1 1 , 

MNP 

♦ 

BVAA I  (20,211  ,VBR(2G,21I  ,Y81  (20, ELI, nrllOOl.YCRI 20,21), V' II 20,211, 

MNP 

5 

CYSRI20) • V  S 1  (201 ,KEEP(20) 

MNP 

t 

REAC  (1,1001  M«(MS(1I,I-1  «M|  «Ncut. 

MNP 

7 

100  FORMAT  (8110  ) 

MNP 

e 

READ  (1,3001  PC TR. PCT BR.PCT  1  iP-.T-ii ,  12 

MNP 

s 

rxITE  (3,1101  PCTR,PCTBR,PCTl,PL(ai,U 

MNP 

1C 

110  FORMAT  Cl*, TIC, • MAX  RANO  ERKOn  uN  REAL  «• F6. 3, 10X, •  B1 AS 

ERROR  MNP 

11 

A  UN  REAL  -*F6. 3»*  OF  ELEMENT*1/ Tlu,*  ON  (MAGINAR  MNP 

12 

BY- <F6.3,21X*0N  IMAGINARY-*F6. 3, is*. *StEO-*  1 5/// 1 

MNP 

13 

REAC  (121  HT .HEAD, NF, NO 

MNP 

14 

N-ND-M 

MNP 

15 

READ  (1,100  1  (  KEEP  1 1 1 , 1 -1,N  1 

MNP 

16 

NERR-C 

MNP 

17 

00  2SC  L-l ,NF 

1MNP 

18 

IX-IZ«2+1 

1MNP 

IS 

REAC  (121  H2(LI,((VR(I  ,JI  ,  VI 1 ! ,01 ,1-1 ,NOI ,j-l, NDI 

1MNP 

20 

(■RITE  (3,120  1  H2(  L) 

1MNP 

21 

120  FORMAT  ( •  1* /T30 , *MA IN  SYSTEM  0,  sPL-ctMENT  MOBILITY  REAL, 

IMAGI  MNP 

22 

A NARY  FREQ-<F8.2//I 

IHNP 

23 

CALL  H0UT2  (VR.NO.NO  1 

1MNP 

24 

CALL  M0UT2  (VI, NO, NO  1 

1MNP 

25 

REAO  (111  HF(  L)  •  (  (  VBR(  I.JI.YBIll.J)  ,i-l,MI,J-l,MI 

1MNP 

26 

WRITE  (3,13i:  HF (L 1 

1MNP 

27 

130  FORMAT  (*1*/T30,'  SUBSYSTEM  01  st'bAcEMENT  MOBILITY  REAL, 

IMAGI  1MNP 

28 

ANAR V  FRE0-*F7.2//I 

IMNP 

29 

CALL  M0UT2  (  VBR.M.M  1 

1MNP 

30 

CALL  M0UT2  (  YBl.M.M  1 

IMNP 

31 

IF ( FE  (L  l-MF (LI)  140,150,140 

IMNP 

32 

140  WRITE  (3,3101 

IMNP 

33 

GO  TC  350 

IMNP 

34 

150  IF(PCTR.NE.O.OK.PCT8R.NE.O.OR.PCli.Nt.O.OR.PCTBI.NE.O  1  NERR- 

1 

IMNP 

35 

IF  (NERR.NE.  1  1  CO  TO  180 

IMNP 

36 

CALL  ERRNU  I V » , VI ,PCTR, PCTBR.Pc 1 1 ,Pd BI ,ND,NO, IX  I 

IMNP 

37 

WRITE  (3,160  1 

IMNP 

38 

160  FORMAT  1* 1*/T15 , ‘MAIN  SYSTEM  OI*PtAceHENT  MOBILITY  WITH  ERROR 

IMNP 

39 

AREAL, IMAGINARY  FREQ-'F8.2//1 

IMNP 

40 

CALL  M0UT2  (  VR.NO.NO  I 

IMNP 

41 

CALL  M0UT2  (  Yl  ,NO,NO  1 

IMNP 

42 

WRITE  13,110  I 

IMNP 

43 

170  FORMAT  ♦,1*/T15,*  SUBSYSTEM  DISPLACEMENT  MOBILITY  WITH  ERROR 

IMNP 

44 

AREAL,  IMAGINARY  FREO*-* F8, 2// 1 

IMNP 

45 

CALL  ERRNU  ( YBR,  YBI  ,PCTR,PCT8M,Pci  1  .PCTBI , NO,  NO,  IX  1 

IMNP 

46 

CALL  M0UT2  (  VBR.M.M  1 

IMNP 

47 

CALL  M0UT2  (  YBl.M.M  ) 

IMNP 

48 

IMNP 

49 

180  00  ISC  1-1, N 

2MNP 

50 

YAARI I, ll-YR(KEEP) 1 1 , N'OL  1 

2MNP 

51 

HO  YAAI 1 1 , ll-Yl ( XEEPI 1 1 ,NCOL  1 

2MNP 

52 

00  2CG  1*1, M 

2MNP 

53 

YCR  1 1 , 1 1  -  YR (MS( 1 1 , NCOL  I 

2  MNP 

54 

200  YCI ( 1 ,1 1 -VI IMS( 1 1 , NCOL  I 

2MNP 

55 

30 


OU  210  l  *  l.N 

2HNP 

56 

OC  21C  J" l » M 

3NNP 

57 

YABRII, JI-VRIXEEPI  1 1« MSI  J)l 

3MNP 

5* 

210 

VABI (l,J)*YI(KE£P(l),HS(J)t 

3MNP 

59 

00  22C  l«l,M 

2MNP 

6C 

00  220  J*1,M 

3MNP 

61 

VCCRI I ,J)*YR(MS(i).MS(JI) 

3MNP 

62 

220 

YCC  l 1 1 1 J 1 *Y1 ( MS ( I  I fMSIJI 1 

3MNP 

63 

IMNP 

64 

230 

K*N*l 

1MNP 

65 

LL«1 

IMNP 

66 

CO  250  I *K , NO 

2MNP 

67 

YAARCI, 11-0. 

2MNP 

68 

YAAHI.ll-O. 

2MNP 

69 

DO  240  J«l,M 

3MNP 

70 

YABP (i,J)*-YBR(LL,J) 

3MNP 

71 

240 

V4BIII, J)*-VB1(LL,JI 

3MNP 

72 

250 

LL-LLU 

2MNP 

73 

CALL  M0LT2  (YAAR.NC.l  i 

IMNP 

74 

CALL  M0UT2  (VAAl.ND.l  1 

IMNP 

75 

CALL  MUUT2  (YABR.NC.M  1 

IMNP 

76 

CALL  M0UT2  (YABI.ND.M  1 

IMNP 

77 

IMNP 

78 

IMNP 

79 

CALL  M0UT2  (VCR.M.l  1 

IMNP 

80 

CALL  M0UT2  (VCI.M.l  1 

IMNP 

81 

CALL  MUUT2  (VCCR-M.M  1 

IMNP 

82 

CALL  HOOT2  (VCCl.M.M  1 

IMNP 

83 

CALL  MATAS  I YBK, YCCK ,M,M, l .  ) 

IMNP 

84 

CALL  MATAS  ( VB1 . VCCI ,MtM,l.  1 

IMNP 

85 

CALL  CINV  1 YBR.YBI .M.YCCR.YCCI  1 

IMNP 

86 

CALL  CMNPY  ( YABR.Y ABI , VCCR.  VCC1 .NO.M.M. VR, VI  1 

IMNP 

87 

CALL  CMMPY  1  YR, Yl « VCR ,YC  I  .NO.M. t, TAttM 

,YA3I  I 

IMNP 

88 

CALL  MATAS  (  YAAR , YABR .NO.l ,-l.  1 

IMNP 

89 

CALL  MATAS  (  YAAI , VABI ,N0 , l ,-1.  1 

IMNP 

90 

WRITE  (3.260)  H2 1 L 1 ,NCCL  ,  (MSI  1 1 «l-i ,M 

1 

IMNP 

91 

260 

FORMAT  I'l'/TlO, 'DISPLACEMENT  MuoUiTV 

REAL,  IMAGINARY  FREQ** 

IMNP 

92 

1 

k  FB.2 .*  HER T2,/T20.* FORCING  AT  MAIN  aYSTEM  ST  A-' I  3/T20, • SUB SV STE M 

IMNP 

93 

BIN 1EHFACL  STA* • 1 3 . • , • , 13, • ,  •  . U//// ) 

IMNP 

94 

WRITE  (3.3201 

IMNP 

95 

00  27C  1*1. N 

2MNP 

96 

270 

WRITE  13.340  )  YAAR  1 1  ,1)  .VAAI(i.i) 

2  MNP 

97 

WRITE  (3.330) 

IMNP 

98 

DC  2EC  1*1, M 

2MNP 

99 

280 

WHITE  1  3.340  )  VAARINaI  •  1 )  ,  YAAl  (Mi  ,13 

2MNP 

100 

290 

CONTINUE 

IMNP 

101 

300 

FURMATI  4F10.4.U0  ) 

MNP 

102 

310 

FORMAT  (•!•//•  FORCING  FREQUEnvIca  OF 

MAIN  BODY  AND  SUBSYSTEM  ARE 

!  MNP 

103 

t  INCOMPATIBLE  •/•  JOB  TEBMINaUo'I 

MNP 

104 

320 

FORMAT  (  20X . 'MAIN  SYSTEM  RESPunac'/) 

MNP 

105 

330 

FORMAT  (  20X,' INTERFACE  RESPJnac'/I 

MNP 

108 

340 

FORMAT  (  20X.1P2E15.4) 

MNP 

107 

350 

CALL  EXIT 

MNP 

108 

END 

MNP 

109 

31 


StttRCUT INfc  MATAS  (  A,B.N1,N2,>  j 
C  AOOIT IOf«  CF  MATRICES  AIN1.N2I  MftO  oiNl*N21  STORED  IN  A 
CIMERSI ON  A 1 20 • l > • B (20 > l > 

DO  100  1-1. Ml 
00  ICC  J-1.M2 
100  AII.JI-AII.JI«S*81I,JI 
RETURN 
END 


MAT  1 
MAT  2 
MAT  1 
1MAT  A 
2MAT  5 
2MAT  A 
MAT  7 
MAT  t 
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SUBROUTINE  ERRNU  (A  ,B  ,  PC  TR,  PUTom  ,Pw  1  i  .PCTBI ,  NJ.NP.IX) 

ERR 

c 

ERR 

c 

A  BIAS  ERROR. 

ERR 

c 

PCTb  (RATIO)  ON  AMPLITjul.  ANO  A  UNIFORM  RANOON  ERROR 

ERR 

c 

HAVING  4  ♦/-  MAXIMUM  0 f  fs.1  (KAMO)  ON  AMPLITUOE. 

ERR 

c 

USES  RANOU 

ERR 

c 

ERR 

DIMEKSICN  AI20.21I .BI20.21) 

ERR 

IFIPCTRI  110.10C.110 

ERR 

100 

IFIPCTBRI  110.130.110 

ERR 

10 

no 

00  12C  1*1. NJ 

1ERR 

11 

00  12C  J*1,NP 

2ERR 

12 

CAU  RANOU  UX.IV.YFLI 

2ERR 

13 

IX*  IV 

2  ERR 

14 

E*l.C«2.0*PCTR»(VFL-O.SI«PCTBR 

2  ERR 

IS 

All  , J)*A 1 1 , J | *E 

2ERR 

1C 

CALL  RANOU  IIX.IV.VFLI 

2ERR 

IT 

IX*IV 

2ERR 

IB 

E-1.C*2.0*PCT1  *(VFL— 0*51 ♦PCI B1 

2ERR 

19 

120 

Bf  I.JI-BI  I.JME 

2ERR 

20 

130 

RETURN 

ERR 

21 

END 

ERR 

22 
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SUBROUTINE  HOUT2  IA,M,NI 
REAL  AI20.21) 

I0-NINOIN,10I 
NR  I  IE  13,1001  (1,1-1,101 
100  FORMAT  1/15,101121 
UNITE  13,1001 
00  UC  1*1, M 

110  URITE  (3,1201  1 • 1 A 1 1 , J ) , J*1 ,1  Cl 
120  FORMAT  1 1 5, 5X,1 P10E 12, 4) 

IF  ( 1D-M  130,170,170 
130  I0*MIN0(N,201 

URITE  (1,1001  (1,1-11,10  I 
URITE  13,1001 
00  1AC  1-1, M 

140  UNITE  13,1201  1 , (A  1 1 , J I , J- 1 1 , 1 «  1 
1FII0-NI  150,170,170 
150  URITE  (3,1001  (1,I-21,N  I 
URITE  (3,100  I 
00  1«C  1-1, N 

160  URITE  13,120  I  1,1  All ,3) ,J-2l  ,ft  I 
170  RETURN 
END 


NOT  1 
MOT  2 
MOT  3 
NOT  A 
MOT  5 
NOT  t 
1M0T  7 
1M0T  8 
MOT  9 
NOT  10 
NOT  II 
NOT  12 
NOT  13 
1N0T  14 
1M0T  15 
NOT  16 
NOT  17 
NOT  18 
1N0T  19 
1M0T  2C 
MOT  21 
NOT  22 
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S  TiSCCT  INE  HNPY  J  A ,  E  ,N1 ,  HZ ,  N3.wi 

NPY 

1 

c 

NPY 

2 

c 

C  •  A  •  a 

NPY 

3 

c 

*  INI  X  ,ui 

ft  (N2  X  Nil 

C  INI  X  N3I 

PRY 

A 

c 

NPV 

t 

REAL  A<20, 211 ,9120,211, CI20.21J 

npy 

6 

00  ICO  I*l,Nl 

1NPY 

7 

00  ICC  J«1,N3 

2NPY 

8 

CII.JI-C. 

2NPY 

< 

00  ICC  X«1,N2 

3NPY 

10 

loo  ciiiJt-c(iiJ)*A(i.K)*a(K,ji 

3NPY 

11 

RETURN 

NPV 

12 

END 

NPV 

13 
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SOBRCLT INfc  INVRS  (B.N.AI 

INV 

1 

A  -  INVERSE  OF  B  B  UNOIjIunotU 

IN- 

2 

INV 

3 

DIMENSION  A  ( 20  >21 )  ,0120, 211, IRu«m!  ,  1  COL  (211,6(20,21 

INV 

4 

00  ICC  I«  1,  N 

1  INV 

5 

DO  ICC  J»1,N 

2I*!V 

6 

1U0 

A(  I >  J)*B( l,JI 

2  INV 

7 

M*N4  1 

INV 

e 

DO  11C  1 *1  > N 

1  INV 

9 

1 KGh  1  11*1 

1  INV 

10 

no 

ICOUIIM 

1  INV 

11 

00  26C  A* 1 >N 

1  INV 

12 

AMAX  *  AIK.KI 

1INV 

13 

OC  1 3C  1  *K ,  N 

2  INV 

14 

00  13C  J*K,N 

3  INV 

1! 

IF1ABS1  A  1 1 > J| l-ABSI AMAX! 1 1 30  > L*<J, 140 

3  INV 

16 

no 

AMAX*  A  1 1 • J) 

3  INV 

1  7 

IC*I 

3  INV 

IB 

JC»J 

3  INV 

19 

130 

CONTINUE 

3  INV 

20 

Kl-ICCLIKI 

1  INV 

21 

ICUL  1  Kl « ICQL ( 1C  1 

1  INV 

22 

iconic  l«Kl 

1  INV 

23 

K  1  *  1  RUM  IK) 

1  INV 

24 

IROtalK) *1RCr( JC) 

1INV 

25 

IKOUIXI-Kl 

1  INV 

26 

1F1APAX!  160.14C.160 

1  INV 

27 

no 

RRITE  I  3.  150 

1  INV 

2e 

150 

FORHATI •  SOLUTION  OF  EXISTING  MmIkia  nUT  POSSIBLE'! 

1  INV 

25 

GO  TC  330 

1  INV 

30 

160 

00  17C  J*1.N 

2  INV 

31 

E-A1K.JI 

2  INV 

32 

AIK, J)*A1 IC.JI 

2  INV 

33 

170 

Al  IC.JI  *E 

2  INV 

34 

CO  1 FC  !«1.N 

2  INV 

35 

t  *  A  1  l.KI 

2  INV 

36 

Al 1  ,Kt*A( I,  JCI 

2  INV 

37 

loo 

Kl  1  .  JCI  *£ 

2  INV 

38 

OC  2  1C  1  *  1 .  N 

2  INV 

39 

I F 1  1-KI  200.190,200 

2  INV 

40 

190 

Al 1 ,F»*1. 

2  INV 

41 

GO  TC  210 

2  INV 

42 

200 

Al 1 ,PI*C. 

2  INV 

43 

210 

CONTINUE 

2  INV 

44 

PVT*A1K,KI 

1  INV 

45 

DO  22C  J*l,H 

2  INV 

46 

220 

AIK, Jl-AIK.JI/PVT 

2  INV 

47 

ou  lie  i*i. n 

2  INV 

48 

IF! 1-KI 2JC.230.230 

2  INV 

45 

230 

arult*au  .x> 

2  INV 

50 

UO  24 C  J*  l » R 

3INV 

M 

24u 

Al I , J)*AI I,JI-ARULI*A1K, Jl 

3  INV 

52 

250 

CONI  INUt 

2  INV 

53 

00  2tC  1*1, N 

2  INV 

54 

260 

A  1  1  .Kl  *A1  1  .Ml 

2  INV 

55 

36 


! 


00  2  So  1*1, N 

1INV 

5* 

00  2 7C  1*1, N 

21NV 

37 

IFI  IkCNlI »-LI270,2»0.270 

2INV 

5f 

270 

CONTINUE 

2 INV 

59 

260 

00  2SC  J* 1 , N 

2INV 

6C 

290 

OIL , J 1 *A{ 1 , J 1 

2  INV 

61 

DO  J2C  J*  1 , N 

1  INV 

62 

00  3CC  L*1,N 

2  INV 

63 

IF  i ICCL 1 J 1-L )  300,310,300 

2  INV 

6A 

300 

CONTINUE 

2  INV 

65 

310 

00  320  I-l.N 

2  INV 

66 

320 

All ,1 l*C 1 1 , J) 

2  INV 

67 

330 

RETOHN 

INV 

6E 

ENO 

INV 

65 

c 


SODROUT  INfc  RANOU  <1*.1Y,YFU 

THIS  SUBROUTINE  a  riiuH  SSP  VERS,  II 

iY«im»3f 
IFI1Y1  ICC. 110, 110 
lUO  I Y» 1 Y*2 147483A4TM 
110  VFL-IV 

YFL»YFL*.<i656tl3E-9 

RETURN 

END 


RAN  1 
RAN  2 
RAN  3 
RAN  A 
RAN  5 
RAN  t 
RAN  7 
RAN  B 
RAN  9 
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guou 


SUBROUTINE  CINV  U.B.N.C.DI 


C*I*0  -  i.W cKit  OF  A»l*B  l-SURTI-U 

B  ASSURcJ  non  SINGULAR 


REAL  AI20.21I  .6 120. 21 1. C <20. 21 1, 0120.211. El  20. 211 
CALL  INVRSI B.N.CI 
CALL  HNPYIC.A.N.N.N.E) 

C*LL  RMPYlA.fc.N.N.N.CI 
00  ICC  1*1. N 
00  ICC  J*1.N 
100  C(I.JI-C(I.J)»B(l.JI 
CALL  INVRSIC. N.CI 
CALL  MPY1E.0.N.N.N.CI 
GO  110  1*1,  N 
DO  110  <J*1.N 
110  OII.JI— OII.JI 
RETURN 
ENO 


1C  IN  1C 
2CIN  11 
2CIN  12 
CIN  13 
CIN  14 
1C  IN  IS 
2CIN  16 
2CIN  17 
CIN  18 
CIN  19 
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c 

c 

c 

c 

c 


SUBRCLT  I.NE  CMMPV  < A,B,C,0,Nt,N*,lw,  t,F  | 


COMPLEX  MATRIX  MULT 
t  ♦  l*F  •II  ♦  l«BI«IC  ♦  l*Oi 
A.B  ARC  hi  X  Mi  CiC  ARE  * 


1  •  SQRI-ll 

HJ  £.F  ARE  N1  X  N3 


DIMENSION 

*.GI2n,21) 


Al 20,21)  #BI2 0t21l  ,CUu,4ii,i,t20,21l,EI20,21l  , FI 20,21) 


CALL  ^MPV  IA,C,Nl«N2,N3,E  ) 
CALL  MNPY  IB,0,N1,N2,N3,G  I 
00  ICC  1-1, N1 
00  ICO  J-I.N3 
100  EII,J)-EII,J)~GII,JI 

CALL  MMPV  I A,0,N1,K2,N3,F  I 
CALL  MMPV  |  8,C,N1,N2,N3,G  ) 
00  110  1-1, N1 
JO  11C  J-l.NJ 
110  FII,J)-FI l,JI*G(l,JI 
RETURN 


ENO 


CMM 

1 

CMM 

2 

CMM 

3 

CMM 

4 

CMM 

5 

CMM 

6 

CMM 

7 

CMM 

6 

CMM 

9 

CMM 

10 

1CMM 

11 

2  CMM 

12 

2CMM 

13 

CMM 

14 

CMM 

15 

1CMM 

16 

2CMM 

17 

2  CMM 

IE 

CMM 

19 

CMM 

20 

* 
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